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We calculate the conductance of a circular constriction of radius a in an insulating diaphragm 
which separates two conducting half-spaces characterized by the mean free path I. Using the Boltz- 
mann equation we obtain an answer for all values of the ratio (/a. Our exact result interpolates 
between the Maxwell conductance in diffusive (I <C a) and the Sharvin conductance in ballistic 
(I S> a) transport regime. Following the earlier approach of Wexler we find the explicit form of the 
Green's function for the linearized Boltzmann operator. The formula for the conductance deviates 
by less than 11% from the naive interpolation formula obtained by adding resistances in the diffusive 
and the ballistic regime. 
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I. INTRODUCTION 

The problem of electron transport through an orifice (also known as a point contact) in an insulating diaphragm 
separating two large conductors (Fig. has been studied for more than a century. MaxwellH found the resistance in 
the diffusive regime when the characteristic dimension a (radius of the orifice) is much larger than the mean free path 
i. Maxwell's answer, obtained from the solution of Poisson equation and Ohm's law, is 



RM = 2? 



(1) 



where p is resistivity of the conductor on each side of the diaphragm. Later on, Sharvinu calculated the resistance in 
the ballistic regime (£ ^> a) 



4p£ /2e 2 kpA 
S 3A { h An 



(2) 



where A is the area of the orifice. 
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FIG. 1. Electron transport through the circular constriction in an insulating diaphragm separating two conducting half-spaces 
(each with a mean free path £). 



This "contact resistance" persists even for the ideal conductors (no scattering) and has a purely geometrical origin, 
because only a finite current can flow through a finite size orifice for a given voltage. In the Landauer-Biittiker 
transmission formalisma, we can think of a reflection when a large number of transverse propagating modes in the 
reservoirs matches a small number of propagating modes in the orifice. In the intermediate regime, when a ~ £, the 
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crossover from Rm to Rs was studied by WexlerB using the Boltzmann equation in a relaxation time approximation. 
The influence of inelastic collisions on the orifice current-voltage characteristics was studied using classical kinetic 
equations in Ref. [s] and quantum kinetic equations (Keldysh formalism) in Ref . |[ This effect underlies an experimental 
technique for the extraction of the phonon density of states from the nonlinear current-voltage characteristics (point 
contact spectroscopy^!) . Recentlj£,_the size of orifice has been shrunk to a ~ Xp allowing the observation of quantum- 
size effects on the conductancdafl. In the case of a tapered orifice on each side of a short constriction between 
reservoirs, discrete transverse states ("quantum channels") below the Fermi energy which can propagate through the 
orifice give rise to a quantum version of Eq. (^). The quantum point contact conductance is equal to an integer 
number of conductance quanta 2e 2 fh. 

Here we report a semiclassical treatment using the Boltzmann equation. Bloch-wave propagation and Fermi- 
Dirac statistics are included, but quantum interference effects are neglected. Electrons are scattered specularly and 
clastically at the diaphragm separating the electrodes made of material with a spherical Fermi surface. Collisions are 
taken into account through the mean free path I. A peculiar feature is that the driving force can change rapidly on 
the length scale of a mean free path around the orifi.ee region. The local current density depends on the driving force 
at all other points. Our approach follows Wexler'aa study. We find an explicit form of the Green's function for the 
integro-differential Boltzmann operator. The Green's function becomes the kernel of an integral equation defined on 
the compact domain of orifice. Solution of this integral equation gives the deviation from the equilibrium distribution 
function on the orifice. Therefore, it defines the current through the orifice and its resistance. 

The exact answer can be written as 



R(e/a)=R s +~f(£/a)R M , 



(3) 



where j(£/a) has the limiting value 1 as l/a — > and Rs/Rm — * 0. We are able to compute j(£/a) numerically to 
an accuracy of better than 1%. Our calculation is shown on Fig. 0. We also find the first order Pade fit 



7fit(Z/ a ) = 



1 + 0.83 l/a 
1 + 1.33 l/a 



which is accurate to about 1%. Our answer for 7 differs little from the approximate answer of Wexler 
on Fig. I as 7 Wcx . 



(4) 

also shown 




FIG. 2. The dependence of factor 7 in Eq. (pOI) on the ratio l/a. Also shown is the variational calculation of Twex from Ref. 



Section || formulates the algebra and Sec. Ill explains the solution 



II. SEMICLASSICAL TRANSPORT THEORY IN THE ORIFICE GEOMETRY 



In order to find the current density j(r) through orifice, in the semiclassical approach, we have to solve simultaneously 
the stationary Boltzmann equation in the presence of an electric field and the Poisson equation for the electric potential 
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. dF(k,r) eV$(r) dF(k,r) _ F(k, r) - / LE (k, r) 

r ' dr h ' 9k r [) 



V ^(r)=-^, (6) 
e 

Mr) = ^£(F(k,r)-/( efc )), ( ? ) 

k 

= ^£(^(k,r)-/ LE (k,r)), (8) 

k 

j( r ) = ^E v ^( k ' r )- ( 9 ) 



Here -F(k, r) is the distribution function, /(e^) is the equilibrium Fermi-Dirac function, <I>(r) is electric potential, £1 
is the volume of the sample and /le (k, r) is a Fermi-Dirac function with spatially varying chemical potential /z(r) 
which has the same local charge density as F(k, r). In general, we have to deal with the local deviation Sn(r) of 
electron density from its equilibrium value self-consistently. The collision integral is written in the standard relaxation 
time approximation with scattering time r = 1/vf- This system of equations should be supplemented with boundary 
conditions on the left electrode (LE) at z = — oo, right electrode (RE) at z = oo, and on the impermeable diaphragm 
(£>) at z = 0: 

$(r LE ) = V, (10a) 
$(r RE ) - -V, (10b) 
3z(r D ) = 0, (10c) 

where z-axis is taken to be perpendicular to the orifice. In linear approximation we can express the distribution func- 
tion F(k, r) and local equilibrium distribution function /le (k, r) using 5{i(r) (local change of the chemical potential) 
and ^(k, r) (deviation function, i.e. energy shift of the altered distribution) 

/le (k, r) = f(e k - ^(r)) « /(e fe ) - ^^(r), (H) 
F(k, r) = f(e k - *(k, r)) » /(e fe ) - ^^*(k, r). (12) 
These equations imply that Sfi(r) is identical to the angular average of ^(k, r) 

Mr) = ^ -^^*( k ' r ) = ^(°)<*( r )) = ^(°)^(r), (13) 

k 

where iV(0) is the density of states at the Fermi energy cf- In the case of a spherical Fermi surface, 

(*(r)) = -?- I dO fe *(k,r). (14) 



4ir _ 

Following Wexleri, we introduce a function u(k, r) by writing ^(k, r) as 

#(k,r) = eVu(k,r) - e$(r). (15) 
Thereby, the linearized Boltzmann equation (^]) becomes an integro-differential equation for the function u(k, r) 

^■^^- = (u(v))-u(k,v). (16) 

To solve this equation we need to know only boundary conditions satisfied by u(k, r) and then we can use this solution 
to find the potential ^(r). Thus the calculation of the conductance from u(k, r) is decoupled from the Poisson equation. 
This is an intrinsic property of linear response theorylld. The boundary conditions for ( |l6| ) are: 

(u(r LE )) = 1, (17a) 
MntE)) = -i. (17b) 
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They follow from the boundary conditions ( |10a| )-( |l0b| ) for the potential $(r) and the fact that far away from the 
orifice we can expect local charge neutrality entailing 

Mr)) - -f 1 - (18) 

The driving force is explicitly absent from (|l6|), but it enters the problem through these boundary conditions. Since 
Eq. ( |l6| ) is invariant under the reflection in the plane of diaphragm 

(k,r)-(k«r*), (19a) 
r R =(x,y,-z), (19b) 
k R = (k x ,k y ,-k z ), (19c) 

the boundary conditions imply that u(k, r) has reflection antisymmetry 

u(k,r) = -u(k R ,r R ). (20) 

Wexler's solutionllto the equation ( |l6| ) relied on the equivalence between the problem of orifice resistance and spreading 
resistance of a disk electrode in the place of orifice. Technically this is achieved by switching from the equation for 
function u(k, r) to the equation for function 

iu(k,r) = l + sgn(z)u(k,r). (21) 

The beauty of this transformation is that new function allows us to replace the discontinuous behavior of u(k, r) on 
the diaphragm (which is the mathematical formulation of specular scattering) 

u(k,r D - v k eft) = u(k R , r D - v k dt) = ~u(k, v D + v k eft) (22) 

with continuous behavior of w(k, r) over the diaphragm, discontinuous behavior over the orifice and simpler boundary 
conditions on the electrodes 

(w(r LE )) = H>re)) = 0. (23) 

The Boltzmann equation ( |l6| ) now becomes 

i k ■ dw( x ,V) + w(k, r) - Mr)) = s(k, r)S(z)9(a - r), (24) 
or 

where we have introduced the function 

s(k,r) = 24 z K(k,r) (25) 
which is confined to the orifice region. It can be related to w(k, r) at the orifice in the following way: 

s(k,r ) = 2|4z|(l -w{k,r -v k dt)). (26) 



It plays the role of a "source of particles" in Eq. (|24|). The notation ro refers to a vector lying on the orifice, that 
is ro = (x,y,0) with x 2 + y 2 < a 2 . The discontinuity of w(k, r) on the orifice is handled by replacing it by the disk 
electrode which spreads particles into a scattering medium. 

The Green's function for Eq. ( p4[ ) is the inverse Boltzmann operator (including boundary conditions) 

4 • ^ + 1 - A\ G B (k, r; k', r') = d(Q k - il k , )S(r - r'), (27) 

and A is the angular average operator 

i/(k) = i- /dO fc /(k)=(/). (28) 

The Green's function for the Boltzmann equation allows us to express w (k, r — v^dt ) in the form of a four-dimensional 
integral equation over the surface of the orifice 
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tu(k,r - v k dt) = /" dn /b /rfr' G B (k,r -v k dt;k' ! r' + v k /^)s(k',r , ). (29) 

The function u>(k, r) is discontinuous over the orifice, so we formulate the equation for this function at points in- 
finitesimally close (dt — > +0) to the orifice. We find the following explicit expression for the Green's function 

G B k,r;k',r' = -\ — — - U ft fe - ft fc > + yH H > . 30 

il^l + iq-£ k \ 47r(l + iq-£ k /) / 

Its form reflects the separable structure of Boltzmann operator, i.e. the sum of operators whose factors act in the 
space of functions of either r or k. However it is nontrivial because the factors acting in k-space do not commute and 
the Boltzmann operator is not normal — it does not have the complete set of eigenvectors and the standard procedure 
for constructing the Green's function from the projectors on these states fails. The first term in ( |30"| ) is singular and 
generates the discontinuity of w(k, r) over the orifice. 



III. THE CONDUCTANCE OF THE ORIFICE 



The conductance of the orifice is defined by 



l = ± = [d^jM 

R 2V 2V K ' 

where the z-component of the current at the surface of the orifice is 

J,(ro)= N{Q J e2V I dn k s(k,r ). (32) 
8ttt J 

The Green's function result ( |30| ) allows us to rewrite Eq. ( p9| ) in the following integral equation for the smooth function 
s(k, r ) over the surface of the orifice 



1 = s(k,r ) 



2|4.| 

where G(k, r ; k', r' ) is non-singular part of the Green's function (|3 



J dn k , dr'o G(k, r ; k', r' )s(k', r' ), (33) 



1 r o£e iq(r o _r ' o) 

G(k,r n ;k',r' n ) = Ida-, ^ -. (34) 

V ' °' ' a > 327r 4 7 H (l + iq-4)(^-arctan^)(l + iq-4') 

The distribution function s(k, r ) has two k-space variables, the polar and azimuthal angles (0 k ,4>k) of the vector 
k on the Fermi surface, and the radius tq and azimuthal angle <po of the point ro on the orifice. Because of the 
cylindrical symmetry, s(k, r ) does not depend separately on <f>k, <j)Q, but only on their difference (f>k — 4>o- This allows 
the expansion 

s(k,r ) = ^a L M(r-o)yLM(0fc,^)e- iM * , (35) 

LM 

and Eq. ( |3^ ) can now be rewritten as 

2^cos0 fc = J2 SL'M'(r o )rL'M'(^,0fc)e" lM ' 0o sgn(cos0 k ) 

L'M' 

+2£ f dn fc , dr'o G(k,r ;k',r' ) cos 0* ^ s VM >(r' )Yv U >(e k ,, <j> k ,)e- iM ' *'° . (36) 

L'M' 

This four dimensional integral equation can be reduced to a system of coupled one dimensional Fredholm integral 
equations of the second kind after it is multiplied by Y^ M (8k, 4>k)e lM ^° and integrated over 9 k , 4>k and 4>q. We also 
use the following identities 
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9i 



.92 



! (L-M + l){L + M + l) 
(2L + l)(2L + 3) : 

1 (L - M)(L + M) 
{2L- 1)(2L + 1) ' 



(37a) 
(37b) 

(37c) 



1 /•Wflfc.fo) 



47r J 1 + iq£ k 



and 



2tt 



e iqr Oe -iM0 o ^ 



271 



3 i9xr„o«i(0o-* 4 ) e -iM*o d(f)0 = 2m M 'J M (qjLro)e- m *< , 



(38) 



(39) 



where q_L is projection of q = q z + q^ in the plane of orifice and Jm{%) is the Bessel function of the first kind. For 
the function /l in (38) we get the following expression 



r°° (—i\- L 1 

JlW ) ( l)' / e-*j L (q£x)dx = —l—Q lA — ). 



(40) 



where Jl(x) is spherical Bessel function and Ql{%) is Legendre function of the second kind. Explicit formulae for 
f L (x) are 



arctan x 



x 

-x + arctan a; 



-3a; + (a; 2 + 3) arctan a; 



4™3 



fa(x 
fi(x 
f 2 (x 
fs(x 
U(x 

OX" 

The final form of the integral equation for Slm(^o) i n the expansion of s(k, r ) is 

iiJ-^S L iSM0 = ^2 c LM,L'M'^MM'SL'M'( r o) + 4 / r„ d?"Q iY LM ,L'M' (r Q , r' )s L < M ' (r'o) 

" T ,'M' T.I Ml J° 



2a; 3 

5a; + (5 + 3a; 2 ) arctan; 
2x~ A 



-^a; 3 - 35a; + (35 + 30.t 2 + 3a; 4 ) arctan x 



(41a) 
(41b) 
(41c) 
(41d) 
(41e) 

(42) 



where the kernel of equation is given by 

K^wM) = l M'-M [ _ l) M + M> r q 2 dg r sin6 dQ qt 2 fL'(qt)Yy M ,(6 q ) jM ^gjj^ ^ 

Jo Jo arctan 

x{l L' + L + l { _ 1) L + l gifL+i{qe)YL+iM{9qH (43) 

Kernel (|43} ) does not depend on <p q so that only 0-part of spherical harmonic Ilm(^) (he. associated Legendre 
polynomial) is integrated. The kernel differs from zero only if L + M has parity different from L' + M' . This follows 
from the fact that the kernel is the expectation value 



K hMl vui(ro,r' ) = (LMM\2t cos0G(k,r o ;kV o )|L'M'M'), 
\LMM) =Y Jj y l {e k ,cj>k)e- imo 

of an odd operator under inversion in the basis of functions \LMM). Their parity is given by 

P\LMM) = (-1) L+M \LMM). 



(44) 
(45) 



(46) 



G 



Exactly under this condition the kernel becomes a real quantity. This means that nonzero slm(^o) are real with the 
property 

slmM = (-1) m sl,-mM, (47) 

ensuring that s(k, r ) is real. The conductance is determined by the (L,M) = (0,0) function soo^o)- The non-zero 
slm(?~o) coupled to it are selected by the condition that L + M is even. This follows from s(k, r ) being even under 
reflection in the plane of orifice. Under this operation, cos 9 k — > — cos9 k , but <j>k, <pQ are unchanged; this means that 
the expansion ( |3q ) contains only terms with L + M even. 

The first term on the right hand side in (13fl) is determined by the matrix element 



clm,l'M' = / de k d(/) k smO k Y{; M (6 k ,(l) k )Y L <M>(Ok,(t>k)sg-n(cos0 k ), (48) 

which is expectation value of sgn(cos#k) in the basis of spherical harmonics. It is different from zero if M = M' 
and L — L' is odd. The states must be of of different parity, as determined by L, because sgn(cos6*k) is odd under 
inversion. 

The system of equations (0) can be solved for all possible ratios of l/a by either discretizing variable tq or by 
expanding SL'M'( r o) m terms of the polynomials in ro 

SLM (ro ) = ^2 fl nLMPn {To ) , (49) 
a 

and performing integrations numerically. The polynomials p n {fo) = X)"=o c i r% are orthogonal with respect to the 
scalar product 



ro dr Pn(r )p m (r ) = S nm . (50) 



o 



The first three polynomials are 



y/2 

fWo) = — , (51a) 



pi(r ) = ; 6r ° 4 =, (51b) 

/ x lQV6(r 2 -fr + ^) 

p 2 (r = — , 51c 
aVlOOa 4 - 288a 3 + 306a 2 - 144a + 27 

Each integral equation in the system (^3) then becomes the matrix equation. Their introduction into inherent matrix 
structure of ( E2h gives the matrix equation for either slm(»"o) a t discretized ro or constants o„lm- Therefore, the 
constants o„lm satisfy the following equation 

UaJ^6 L1 6 MO S n0 = ^2 c lm,l'm anL'M + 4 if n ° L L M M <Vi/M', (52a) 

" V n'L'M' 

t^ii'L'M' _-M'-Mf t\M+M' f°° 2 7 C ■ a r in I? 2 fL'jqtyYh'M' (0q) .„ / . „ x .„' , . „ s 
K nhM — l (-1) J 1 d <l J smUqdVq — g £ _ arc t a n q£ — jM ^ q Sm(> q)3M'{ ( l a SVO-Uq) 

x( J L ' +i +H-l) i+ V/l + lM>i+l,Mft)+l i ' + ^ 1 (-l) i " 1 92/L-l(?^-l,M(^)), (52b) 

j^jiqa sm9 q ) = / r dr p n (r )J M (qr sin q ), (52c) 
Jo 

which simplifies using the following result 



.„ . . fl , A a 2+M+ '(g smg^%FKl + f + i; 2 + f + i, 1 + M; -|(ga sin^) 2 ) 

Hfiqa sm&a) — > cj , - , A tt-*- — - , (53) 

ql ^ 2 1 + M (l + f + |)r(l + M) ' V ; 



i=0 



where li^a; 01, 02', z) is a hypergeometric function. The lowest order approximation for s(k, r ) is obtained by 
truncating the polynomial expansion to zeroth order (i.e. constant — which is the space dependence of Sharvin limit) 
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and expansion in spherical harmonics to L = 0. Then the conductance is determined only by the constant aooo 
following trivially from (|52] ) 

N(0)£e 2 a 2 n 

where the lowest order part of the kernel i^oio depends on £/a, 

^nnn 4£ f°° , f ,„ arctangi ( -Zqt + (q 2 £ 2 + 3) arctang^ . „ 2 „ , arctangA (Mqa sm9„)) 2 
KZ= *L ""Jo %l-arctl g H 4?^ (l-3cos^) + — ^ sing, ' 

(55) 

Further corrections are obtained by solving the matrix equation (p3) where the infinite matrix is approximated by 



its finite block. The matrix elements if" L y M (52t) are tedious to compute, but the conductance converges rapidly 



for large n and L. On the other hand, the matrix elements clm,L'M' ( ff8|) are easy to compute but the conductance 
converges slowly in the ballistic limit which is determined by this matrix elements. We keep the low order matrix 
elements ]f° L y M but go to high order in Clm.i/m'- In practice we find that for c-matrix L max = 12 is sufficient, 
whereas for if-matrix the approximation L max = 2, n max = 2 gives convergence to 1%. The conductance as a 
function of £/a is shown on Fig. ^|. It is normalized to the Sharvin conductance, i.e. conductance in the limit £ 3> a, 
for which 

G(k,r;k',r')-0, 

s(k,r) = 2|4 Z |- (56) 



In the opposite (Maxwell) limit, when £ <C a, we have 

q£ 3 



q£ — arctan ql (q£) 



2 



9/5 + o((q£Y), (57a) 



G(k ' r;k ' r ' } / dq£ W^ = lSir^r-rT ^ 



which is the standard Green's function for the Poisson equation. The dependence of the full Green's function (|30|) on 
k vector is reflection of non-locality. The conductance in the transition region from Maxwell to Sharvin limit can be 
compared with the naive interpolation formula which approximates resistance of orifice by the sum of Sharvin and 
Maxwell resistances 

<5-*-*( 1+ t!)- < 58 ' 

Somewhat unexpectedly, naive interpolation formula Gj deviates from our result for G at most 11% when £/a — > 1 
as shown on Fig. [|. We can also cast our lowest order approximation for the conductance (|54| ) in an analogous form 
as © 

1 _ /3 32 37ra\ 
7?- 7 + ?Z57-5-7 • (59) 



do V 4 3tt 2 ' 8 £, 

The numerical coefficients in Eq. ( |59| ) are not accurate in this simplest approximation. Replacement of 3/4 by 1 and 
32/(37r 2 ) by 1 yields correct limiting values of G and leads to a plausible interpolation formula. It differs from Eq. ( |5q ) 
by the introduction of a factor 7 which multiplies the Maxwell resistance 



G 



1 „ / 37ra\ 

R s U + 7^7 . ( 6 °) 







tx£ 
16a 



KZ- (61) 



This formula is compared to G and Gi on Fig. It differs from our most accurate calculation of G by less then 1%. 
Therefore, for all practical purposes it can be used as an exact expression for the conductance in this geometry, and 
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it is the main outcome of our wprk. The factor 7 is of order one and depends on the ratio lj a as shown on Fig. 
We also plot on Fig. || Wexler'su previous variational calculation, 7wcx- 




1/ a 

FIG. 3. The conductance G (L = 2, n = 2), normalized by the Sharvin conductance Gs (^), plotted against the ratio i/a. 
It is compared to the naive interpolation formula Gi (^), and the plausible interpolation formula Go (po|). 

In conclusion, we have calculated the conductance of the orifice in all transport regimes, from the diffusive to 
the ballistic. The altered version ( |60| ) of the simplest approximate solution of our theory ([34]) is already reasonably 
accurate. It shows the microscopic theory correction to the naive interpolation formula (sum of Maxwell and Sharvin 
resistances) which is within 11% different from the exact answer. Further corrections converge rapidly to an exact 
result. This analysis is of interest in any situation where the geometry of the sample can introduce additional 
resistances, like for example in transport in granular materialsliil. 

This work was supported in part by the NSF grant no. DMR 9725037. 
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